Gauss's Mathematics
Gauss's work is everywhere in modern physics. The Gaussian integral is the only integral we can do exactly in infinite dimensions—making quantum field theory possible.
The normalized Gaussian defines probability and quantum uncertainty:
Gauss's Law for electricity connects flux to charge—later a constraint in quantum gauge theory:
For surfaces, Gaussian curvature $K$ and the Gauss-Bonnet theorem link local geometry to global topology—central to string worldsheets:
Interactive Demo 1: The Gaussian
Vary $a$ in $e^{-ax^2}$. The area is always $\sqrt{\pi/a}$—the cornerstone of all Gaussian path integrals.
This exact result lets Feynman and Hawking evaluate $\int \mathcal{D}x\, e^{iS/\hbar}$ for quadratic actions.
Maxwell's Theory
Maxwell unified electricity and magnetism. His equations imply light is an electromagnetic wave with $c = 1/\sqrt{\mu_0\varepsilon_0}$.
Modern form uses the gauge potential $A_\mu = (\phi/c, \mathbf{A})$ and field tensor:
Gauge freedom $A_\mu \to A_\mu + \partial_\mu \Lambda$ leaves $F_{\mu\nu}$ invariant. This $U(1)$ symmetry is the template for the Standard Model and for string theory's Chan-Paton factors.
Following Euler’s variational methods and Hilbert’s operator formalism, Maxwell’s theory becomes the first quantum field theory upon quantization.
Interactive Demo 2: Gauss's Law Flux
Flux through any sphere around a charge is $Q/\varepsilon_0$, independent of radius.
In Quantum Mechanics
Quantum states live in Hilbert space. The Gaussian wave packet saturates Heisenberg:
Feynman's path integral is defined by Gaussian integrals. For a free particle:
All evaluations reduce to $\int e^{-a x^2 + bx}dx$. Gauss makes quantum mechanics computable.
Gauss's law as constraint: Physical states satisfy
removing unphysical polarizations—exactly as Dirac constrained Hilbert spaces.
Interactive Demo 3: Wave Packet Spreading
Wider initial packets spread slower. This is quantum diffusion.
In String Theory
String theory’s Polyakov path integral is Gaussian in the embedding fields $X^\mu$:
Evaluating it gives the propagator $\langle X(z)X(w)\rangle \sim -\alpha'\ln|z-w|$, again a Gaussian integral.
The action contains the Gauss-Bonnet term:
Thus amplitudes sum over genus $g$ with weight $g_s^{2g-2}$—topology controls quantum corrections.
Maxwell from strings: Open string massless mode is $A_\mu$. Its low-energy action is Born-Infeld, reducing to
D-branes source RR fields, obeying Gauss's law: $\oint_{S^{8-p}} *F_{p+2}=Q_{\text{Dp}}$.
Interactive Demo 4: Worldsheet Genus
Gauss-Bonnet: $\chi = 2-2g$ determines string coupling weight.
Connections
Gauss → Path Integrals
The Gaussian integral → finite-dimensional Gaussian integrals → infinite-dimensional path integrals. Hawking's Euclidean quantum gravity uses the same $\int e^{-S_E}$.
Maxwell → Gauge Theory
$U(1)$ gauge symmetry generalizes to $SU(3)\times SU(2)\times U(1)$ in the Standard Model, and to non-abelian gauge fields on D-branes in string theory.
Gauss-Bonnet → Topology
Gauss's curvature theorem controls string loop expansion. Euler characteristic, refined by Riemann and Chern, classifies Calabi-Yau compactifications.
GNU Octave Laboratory
Run these in Octave to verify the mathematics Gauss and Maxwell built. All scripts are self-contained and use only core functions.
a) Gaussian Integral Verification
Compare analytic $\sqrt{\pi/a}$ with numerical quadrature.
% Gaussian integral
a = 2;
I_analytic = sqrt(pi/a);
I_numeric = integral(@(x) exp(-a*x.^2), -Inf, Inf);
fprintf('Analytic: %.10f, Numeric: %.10f, Error: %.2e\n', ...
I_analytic, I_numeric, abs(I_analytic-I_numeric));
x = linspace(-3,3,400);
plot(x, exp(-a*x.^2), 'LineWidth',2); grid on;
title('Gaussian e^{-ax^2}'); xlabel('x');
b) Gaussian Wave Packet Evolution
Initial minimum-uncertainty packet $|\psi|^2$.
% Wave packet at t=0
x = linspace(-5,5,500);
sigma0 = 0.5;
psi = (2*pi*sigma0^2)^(-0.25) * exp(-x.^2/(4*sigma0^2));
prob = abs(psi).^2;
plot(x, prob, 'LineWidth',2); grid on;
xlabel('x'); ylabel('|\psi|^2');
title(sprintf('Gaussian packet, \\sigma_0 = %.2f', sigma0));
% Uncertainty: Delta_x * Delta_p = hbar/2
c) Gauss's Law Flux
Flux through sphere independent of radius.
% Gauss's law
Q = 1.602e-19; % elementary charge
eps0 = 8.854187817e-12;
r = logspace(-9,-7,5); % various radii
E = Q ./ (4*pi*eps0*r.^2);
flux = 4*pi*r.^2 .* E; % = Q/eps0
fprintf('Q/eps0 = %.3e\n', Q/eps0);
disp([r' flux']); % constant column
d) Maxwell Wave
Plane wave solution $E = E_0\sin(kx-\omega t)$.
% EM wave
c = 299792458; E0 = 1;
k = 2*pi; omega = c*k;
x = linspace(0,2,500);
t = 0;
E = E0*sin(k*x - omega*t);
B = E0/c * sin(k*x - omega*t);
plot(x,E,'r',x,B*3e8,'b','LineWidth',1.5); grid on;
legend('E (V/m)','cB (V/m)'); xlabel('x (m)');
title('Maxwell wave satisfies \nabla^2E = \mu_0\epsilon_0 \partial_t^2E');
e) Gauss-Bonnet on Sphere
Verify $\int K dA = 4\pi$ for $S^2$.
% Gauss-Bonnet sphere
R = 1;
K = 1/R^2; % Gaussian curvature
area = 4*pi*R^2;
integral_K = K * area;
chi = integral_K/(2*pi); % Euler characteristic
fprintf('∫ K dA = %.5f = 2πχ, χ = %.1f\n', integral_K, chi);
% For genus g: chi = 2-2g, torus g=1 gives 0